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摘 要 影响 热电 制冷 器 仿真 温度 的 输入 参数 较 多 且 不 易 准 确 测 得 导致 建立 高 精度 仿真 模型 较为 困难 . 为 此 , 提出 了 一 
种 热电 制冷 器 热 模 型 建立 与 修正 方法 , 首先 建立 热电 制冷 器 温度 场 仿真 模型 ， 并 基于 Spearman 相关 级 数 对 影响 仿真 温度 
的 参数 进行 灵敏 度 辨 识 ,确定 强 相 关 参 数 和 弱 相 关 参 数 ; 然 后 以 热电 制冷 铝 的 实测 温度 值 为 目标 ,采用 多 岛 群 遗 传 算法 对 
强 相关 参数 进行 修正 ; 在 此 基础 上 ,再 对 弱 相 关 参 数 进 行 修正 ; 最 后 运用 霍 克 基 维 斯 直接 搜索 算法 对 修正 后 的 参数 进行 更 
为 精确 的 修正 ， 得 出 了 电流 为 3.6 A 时 模型 输入 参数 的 修正 值 . 采用 相同 方法 获得 了 电流 在 0~3.9 A 范围 内 输入 参数 的 
修正 值 . 据 此 建立 的 热 模型 ， 仿 真 值 与 实测 值 最 大 误差 小 于 4.9 "C (-5 ~ 100.5°C 范围 内 )， 表 明 该 热 模 型 的 建立 与 修正 
方法 合理 有 效 。 
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Abstract As the input paramete and not easy to measure accurately, it is difficult to get 


the high-precision simulation thermoelectric cooler. A method of model establishment and 
modification was proposed, 位 st， the thermal model of thermoelectric cooler was established. 
Then, the sensitivity of jkermial model Parameters is analyzed by the Spearman rank correlation. It 
is aimed to determine Ne strong correlation parameters and weak correlation Parameters that affect 
the simulation accuracy. Taking the measured temperatures as the fargets, the multi-island genetic 
algorithm was used to correct the strong correlation paramete@ stly. Then the weak correlation 
parameters are corrected. Based on above results, the Hoo -Jeeves Direct Search Method was used 
to modify the parameters for obtaining the accurate va hen the current is 3.6 A. Similarly, the 
‘ved from 0 A to 3.9 A. The accurate thermal 
model of thermoelectric cooler can be established\at different current. The results show that the 
maximum error between simulation value Xi periment value is less than 4.9 °C (—5°C to 100.5 
°C) when the current is changed from 0 XE 3.9 A. In conclusion, the modified method of thermal 


model is reasonable and accurate. 


corrected values are obtained when the current i 
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热电 制冷 右 (TEC) 是 一 种 利用 半导体 材料 热电 
效应 实现 加 热 、 降 温 、 控 温 等 目的 的 器 件 ， 其 广泛 
应 用 于 差 热 发 电 、 半 导体 激光 器 制冷 、 空 调制 冷 等 
领域 . 目前 对 于 热电 制冷 器 的 研究 主要 集中 在 如 何 
有 效 提高 热电 制冷 器 的 优质 系数 [1,a 及 电热 转换 效 
率 8-9. 在 热机 设计 过 程 中 , 建立 准确 的 可 预测 温度 
场 分 布 的 有 限 元 模型 (以 下 简称 “ 热 模型 ”) 是 实现 热 
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电 制 冷 噩 高 效应 用 及 优化 设计 的 基础 热电 制冷 器 
温度 场 的 仿真 本 质 是 热 场 与 电场 的 夭 合 分 析 ， 具 有 
高 度 非 线性 的 特点 ， 获 得 精确 解 较为 困难 . 文献 [7] 
根据 热电 制冷 器 的 工作 原理 结合 实际 工作 情况 ， 推 
导出 了 半导体 制冷 器 在 任意 时 刻 表 面 温度 的 计算 公 
式 , 热 端 实测 值 与 计算 值 最 大 误差 为 12°C. 文献 [8] 
在 考虑 各 种 影响 因素 所 造成 的 附加 热 阻 和 接触 热 阻 
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的 基础 上 , 建立 了 热电 制冷 器 数学 模型 , 获得 了 等 效 
电路 模型 , 并 对 其 特性 进行 了 理论 分 析 . 文献 [9] 建 
芯 了 不 同 散热 方式 下 热电 制冷 器 散热 模型 ， 并 进行 
了 实验 验证 , 最 大 误差 达到 14 "C. 文献 [10], [11] 通 
过 实验 , 研究 了 输入 参数 对 热电 制冷 器 工作 性 能 的 
影响 , 为 后 续 建 立 热电 制冷 妖精 确 模 型 提供 了 思路 。 
现 有 研究 工作 已 为 热电 制冷 右 的 高 效应 用 及 优化 设 
计 做 出 卓有成效 的 贡献 ， 但 模型 的 精度 仍 需 进一步 
提高 . 文献 [12] 运用 混合 算法 对 卫星 的 热 模型 进行 
了 修正 , 修正 后 的 模型 精度 有 所 提高 , 但 是 并 未 将 该 
方法 应 用 于 热电 制冷 器 热 模型 修正 中 . 文献 [13] 建 
立 了 半导体 制冷 系统 中 电 侦 极 在 电场 与 温度 场 的 克 
合作 用 下 非 稳 态 传 热 微 分 方程 ， 通 过 实验 验证 ， 误 
差 到 达 10 "C.。 文献 14 通过 联合 仿真 的 方法 修正 了 
发 热电 阻 散热 模块 ， 修 正之 后 最 大 误差 为 5.4%。 鉴 
于 此 ， 本 文 开展 了 热电 制冷 器 热 模型 的 建立 与 修正 
技术 的 研究 工作 . 

本 文 以 9501/071/040 型 热电 制冷 器 为 研究 对 


象 ， 提出 了 一 种 有 效 提高 热 温度 场 仿真 精度 的 模 殉 六 \ 


在 
能 影响 热 模型 仿真 精度 的 参数 进行 录 敏 及 
es he 

数 作为 决定 模型 仿真 精度 的 输入 和 运用 全 局 优化 算 
法 (Multi-Island Genetic Algofithm) 和 直接 寻 优 算 
法 (Hooke-Jeeves Direct Search Method)， 并 结合 实 
验 数据 对 修正 参数 进行 修正 . 通过 多 次 修正 , 得 到 简 
化 后 仿真 模型 的 输入 参数 关于 电流 的 变化 曲线 ， 最 


1 热电 制冷 右 热 仿真 模型 建立 与 修正 

技术 

热 模型 建立 与 修正 技术 是 将 热 温度 场 仿真 和 测 
试 同时 进行 ， 基 于 修正 算法 ， 使 仿真 温度 与 实测 温 
度 趋 于 一 致 ， 最 终 建立 可 以 精确 预测 不 同 工 况 下 温 
度 场 的 热 模型 . 其 原理 如 图 1 所 示 . 
1.1 热电 制冷 器 误差 来 源 分 析 

热电 制冷 器 是 将 不 同 热电 材料 首尾 串联 ， 利 用 
热电 效应 实现 热量 转移 的 控 温 器 件 . 如 图 2 所 示 为 
热电 制冷 器 的 结构 ， 


吸 热 ( 冷 面 ) 


Cy 


放 热 ( 热 面 ) 
图 2 热电 制冷 右 结 构 示意 图 


Fig. 2 Schematic diagram of thermoelectric cooler 


半导体 材料 的 热电 效应 包含 赛 贝 克 效 应 、 帕 尔 
贴 效应 、 汤姆 还 立 、 焦 耳 效 应 、 健 里 叶 效 应 等 , 表 
达 式 如 (QR. 


发 根据 该 曲线 建立 了 准确 的 热 模型， 全 ja 
RE Qi1=7:T 
Qs =7.71. AT (1) 
AN Q3=1:.R 
BS ND us 


度 辨 识 


图 1 热 模型 建立 与 修正 方法 示意 图 


Fig. 1 Schematic diagram of thermal modeling and 


modification 


其 中 : 马 为 赛 贝克 电 动 势 ; a 为 赛 贝 克 系 数 中 ; Qi 
为 帕 尔 贴 热 ; 了 为 电流 值 ; 7 为 由 尔 贴 系数 ; 82 为 汤 
姆 示 热 ; 7 为 汤姆 逊 系 数 ; Qs 为 焦耳 热 ; Q4 为 傅 里 
叶 热 . 

赛 贝 克 电 动 势 已 与 赛 贝克 系数 a 成 正比 , 帕 尔 
贴 热 Qi 取决 于 电流 值 了 与 帕 尔 贴 系数 rx, 汤姆 进 
热 Q2 与 结 点 温差 、 电 流 值 了 及 汤姆 逊 系数 7 有 关 ， 
这 三 者 代表 了 在 热电 制冷 器 内 部 电 、 热 的 转换 关系 。 
焦耳 热 Qs 可 使 热电 制冷 器 的 冷 热 端 产生 升温 , 傅 里 
叶 热 Q4 决定 了 热电 制冷 器 冷 热 端的 热量 交换 程度 ， 
此 外 ,热电 制冷 器 与 环境 的 热 对 流 以 及 热 辐 射 也 对 
其 温度 场 分 布 产 生 影响 ， 以 上 环节 共同 决定 了 热电 
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制冷 器 温度 分 布 。 

由 此 可 知 ， 热 电 制 冷 器 温度 场 仿真 计算 的 精度 
取决 于 上 述 环节 中 输入 的 参数 是 否 准确 。 半导体 材 
料 物性 参数 随 生 产 批 次 、 工 艺 不 同 、 工 作 温度 变化 
而 不 同 ; 仿真 模型 的 边界 条 件 、 初 始 条 件 、 热 载荷 
等 输入 参数 较 其 真 值 有 所 差异 ; 模型 存在 简化 误差 ; 
以 上 这 些 因素 综合 导致 仿真 结果 与 实测 结果 误差 较 
大 . 热 模型 建立 与 修正 的 原理 就 是 通过 优化 算法 对 
热 仿真 结果 影响 较 大 的 输入 参数 (不 同 模 型 待 修正 
输入 参数 不 同 ) 进行 修正 , 最 终 建立 一 个 能 准确 预测 
仿真 温度 场 的 有 限 元 模型 ， 

1.2 热 模型 参数 灵敏 度 辨识 

由 上 可 知 影响 热 仿真 精度 的 参数 较 多 ， 且 影响 
程度 不 同 。 因 此 需要 在 修正 之 前 进行 参数 灵敏 度 分 
析 ， 确 定 主要 待 修正 参数 . 在 对 比 不 同 灵敏 度 分 析 
方法 的 修正 基础 上 ， 本 文采 用 了 基于 拉丁 超 利 群 抽 
样 的 Spearman 相关 性 灵敏 度 分 析 法 . ee 
关 级 数 可 用 来 考察 不 同 输入 参数 (zizazaz47z5 .， 


响应 量 (温度 值 ) 之 间 的 线性 关系 ， i 


如 式 (2) 所 示 . 


D3 Rr — Q7) o> 
(2) 


rs 二 1 一 
SG 13 一 ty 


式 中 :代表 和 要 ，Ri 9 代表 参数 在 次 
抽样 中 ,第 ;次 样本 点 的 输入 排序 号 、 响 应 排序 号 ， 
当 |rs| > 0.8 时 ， 称 为 高 度 相 关 ; 当 0.8 > |rs| > 0.3 
浆 为 强 相关 ; 当 |rs| < 0.3 称 为 弱 相 关 ， 
1.3 热 仿真 模 型 修正 算法 

热电 制冷 器 的 热 仿真 具有 多 峰值 ， 高 度 非 线 人 
的 特点 . 目前 常用 的 热 模型 修正 算法 有 全 局 优化 法 


名人 化 公法 和 于 等 其中， 全 局 人 人 


法 (如 多 岛 群 遗传 算法 ) 修正 效率 较 低 , 在 迭 

较 低 时 修正 精度 较 差 , 有 旦 提高 有 限 . 梯度 优化 算法 修 
正 效果 较 差 、 易 出 错 , 且 容 易 陷入 局 部 收敛 , 导致 精 
度 不 高 . 直接 搜索 算法 易 陷入 局 部 收敛 , 精度 很 大 程 
度 上 取决 于 初始 值 。 为 提高 模型 修正 的 精度 以 及 收 
敛 特性 ， 本 文采 用 了 多 岛 群 遗 传 算法 进行 全 局 寻 优 ， 
在 得 到 全 局 最 优 解 的 基础 上 ， 以 全 局 最 优 解 为 初始 
值 , 利用 Hooke-Jeeves 直接 搜索 算法 进行 二 次 修正 ， 
既 克 服 了 全 局 优化 算法 效率 低 的 问题 ， 又 提高 了 修 
正 热 模型 的 精度 ， 


2 热 模型 修正 及 验证 


2.1 研究 对 象 建 模 
热电 制冷 器 选用 Ferrotec 9501/071/040 型 ， 最 
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大 许 用 电压 9.8 V. 由 于 热电 制冷 器 在 热 端 散热 较 差 
时 ， 易 造成 温度 过 高 烧 坏 器 件 。 为 此 在 热 面 安装 铝 
制 散 热 器 与 轴 流 式 风扇 ， 冷 面 加 装 小 型 散热 器 ， 强 
化 对 流 换 热 ， 避 免 器 件 烧 坏 .。 实验 对 象 的 几何 模型 
如 图 3(a) 所 示 。 


硅 脂 


人 
热 端 


散热 需 


风扇 


(a) 几何 模型 


(b) 有 限 元 模型 
图 3 热电 制冷 铝 散 热 结 构 模 型 


Fig. 3 Thermoelectric cooler heat dissipation structure model 


© 

建立 有 限 元 模型 ,半导体 BizTes 采 

lectric226 单元 , 铜 导体 、 陶 瓷 片 、 导 
、 铝 制 散热 右 均 采用 热 单元 solid90。 由 于 
作 温 度 较 低 (小 于 100 "“C)， 热 辐射 量 较 对 流 


用 AN 
有 Tm 


NS 故 将 其 忽略 .采用 导热 硅胶 片 连接 各 接 


触 面 强化 热传导 ， 忽 略 其 接触 热 阻 等 。 有 限 元 仿真 
结果 如 图 3(b) 所 示 . 取 热 端 散 热 器 中 间 肋 片上 点 a、 
热 端 与 导热 硅胶 片 共 面 点 5、 冷 端 与 导热 硅胶 片 共 面 
点 c、 冷 端 散热 器 中 间 肋 片上 一 点 4 作为 测 温 点 , 如 
图 4 所 示 ， 

2.2 试验 建立 及 测试 

根据 散热 结构 模型 设计 试验 对 象 ， 搭建 试 验 环 
境 . 温度 采集 及 测量 系统 如 图 4 所 示 . 

温度 采集 及 测 基 系 统 包 括 高 性 能 计算 机 、M+P 
可 开发 动态 信号 识别 分 析 系 统 (温度 采集 模块 )、T 型 
热电 侦 。 

M+P 温度 采集 模块 可 以 同时 测量 48 个 点 的 温 
度 , 采样 精度 达到 0.01 °C, 温度 采集 范围 为 -100 ~ 
十 1400°C, 采样 频率 1000 Hz. T 型 温度 传感器 为 未 
经 封装 的 进口 感 温 线 ， 测 温 精 度 0.1 "C， 测 温 范围 
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T 型 热电 信 电源 由 
图 4 测 温 点 的 布置 以 及 测 温 过 程 示意 图 
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Fig. 4 Experimental installation 


为 -50 ~ 十 200 °C, 通过 对 比 常见 的 2 mm K 型 热电 
偶发 现 , 直径 0.08 mm 的 工 型 热电 偶 裸 线 用 于 测 温 
可 大 大 降低 其 测 温 误 差 .放置 于 热电 制冷 器 结构 中 
可 最 大 程度 减 小 对 热量 传递 的 影响 ， 确 保温 度 采集 
迅速 、 精 确 . 通过 导热 硅胶 片 固定 热电 偶 , 既 实现 传 
感 器 与 被 测 表面 直接 接触 ， 也 降低 了 采用 烙 结 剂 等 
带 来 的 测 温 误 差 . 

将 测试 对 象 放置 于 泡沫 绝热 箱 内 ， 箱 体 四 周 及 
底部 全 部 密封 , 顶部 以 薄板 覆盖 , 降低 外 界 环流 气流 \\ 
re 
架 支 撑 以 降低 热传导 ; 热电 制冷 器 通 直 流 流 
式 风扇 加 载 9.8 V 额定 电压 ; 在 与 仿 相对 应 
的 4 个 点 处 布置 热电 偶 ， 同 时 在 济 浅 寻 架 附近 悬空 
布置 热电 偶 , 测量 环境 温度 . 中 ， 若 5 min 
之 内 温度 值 变化 小 于 0.5*Q\ 助 而 视 为 达到 稳 态 ,以 
此 时 温度 值 为 实验 温度 值 。 利用 实验 测 得 的 温度 值 
对 仿真 模型 进行 修正 及 验证 . 

2.3 参数 灵敏 度 分 析 

影响 热 仿真 精度 的 输入 参数 较 多 ， 需要 对 其 进 
行 灵敏 度 辨识 ,运用 修正 算法 对 灵敏 度 较 大 的 输入 
参数 进行 修正 . 灵敏度 分 析 是 在 输入 参数 的 合理 
值 区 间 内 寻找 最 优 值 ， 需 事先 给 定 所 有 输入 参数 和 
见 的 取 值 范围 . 根据 已 有 研究 结果 15-171, 模型 中 各 
参数 的 取 值 范围 如 表 1 所 示 。 

对 流 换 热 是 被 动 散热 的 主要 散热 方式 ， 对 于 温 
度 场 的 分 布 影响 很 大 . 设 热 端 散 热 器 、 冷 端 散热 器 的 
对 流 换 热 系数 为 G1、G2, 根据 一 个 大 气压 下 、 温度 
为 17 °C 时 的 空气 特性 参数 ， 经 计算 可 知 热 端 换 热 
器 的 对 流 换 热 系数 为 31.4 W/(m? .°C), 于 是 可 将 其 
取 值 范围 设置 为 25~60 W/(m?.°0). 同 理 , 冷 端 对 流 


换 热 系数 的 取 值 范围 设置 为 5~10 W/(m? "CO)D 
表 1 输入 参数 取 值 范围 


Table 1 Uncertain parameters of material 


电阻 率 1/ 导热 系数 2/ 赛 贝克 系数 3/ 


人 10-60.m Wm-1.K-! (10-4V/K) 


(A1):2~6 (A2):0.1~1.5 (As):-4~-1.2 
(B1):2~6 (B2):0.1~2 (B3):0.7~3 
SS (C1):1~2 (C2):350~400 
陶瓷 D (D):15~30 
导热 硅胶 片 媚 (E):1~3 
铝 材 下 (F):150~240 


灵敏 度 分 析 的 第 一 步 是 产生 样本 点 ,样本 点 取 
值 范围 如 表 1 所 示 . 采用 拉丁 超 利 群 法 生成 样本 点 ， 
共生 成 样本 点 300 组 . 将 样本 点 从 阵 加 载 到 模型 中 
进行 计算 , 得 到 仿真 结果 。 假定 4 个 测 温 点 的 温度 
值 为 X。、Xs、XECXKXa, 与 仿真 模型 相同 位 置 处 实验 
对 象 4 个 测 点 的 实验 值 为 ze、zp、zec、za。 修正 目标 
函数 为 实验 得 与 仿真 值 的 最 大 差 值 , 如 式 (3) 所 示 . 
Object=max abs (Xi; — zi) (3) 


2 
\ 


S$。 其 中 , i 为 a, 5, c,d。 编制 Matlab 灵敏 度 算法 , 得 


到 目标 函数 与 所 有 输入 参数 之 间 的 灵敏 度 值 。 根据 
1.2 节 分 析 结 果 将 灵敏 度 进行 排序 , 如 表 2 所 示 。 

表 2 中 , 每 一 项 元 素 与 表 1 中 的 参数 一 一 对 应 ， 
如 : 42 表示 N 型 半导体 导热 系数 . 从 表 2 可 以 看 
出 ， 相 关 性 较 大 的 参数 为 N 型 半导体 赛 贝 克 系数 、 
导热 硅胶 片 导热 系数 、P 型 半导体 赛 贝 克 系 数 ,P 型 
半导体 导热 系数 、 热 端 对 流 换 热 系数 、N 型 半导体 导 


表 2 参数 灵敏 度 结果 


Table 2 Parameter sensitivity analysis result 


所 有 参数 按 灵敏 度 从 大 到 小 排序 


As 五 万 3 万 2 G1 42 G2 
一 0.48 一 0.34 0.33 0.25 0.25 0.13 0.12 


Bi1 D F C1 41 C2 
一 0.12 一 0.09 0.07 0.06 0.05 0.02 


Chin axX1 人 V- 
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因此 ， 可 设置 A43、 忆 、Bs 为 强 相关 参数 ， 到 的 参数 作为 初始 值 , 采用 Hook-Jeeves 直接 搜索 算 


Ba。、G1、A2、G2 为 弱 关 参数 ， 其 他 参数 不 作为 修 
正 参 数 ， 
2.4 模型 修正 及 实验 结果 分 析 

设置 热 模型 的 仿真 条 件 与 实验 条 件 一 致 ， 计 算 
工作 电流 为 3.6 A 时 ,修正 前 热电 制冷 器 4 个 测 点 
温度 值 , 并 与 实验 值 进行 对 比 , 如 表 3 所 示 ， 


表 3 3.6 A 时 测 点 温度 仿真 值 与 实验 值 的 对 比 


Table 3 Simulation value and experimental value 


comparison 
测 点 a/°C b/°C c/°C d/°C 
仿真 值 28.69 69.13 23.92 25.0 
实验 值 36.5 84.6 3.8 5.5 
误差 一 7.81 15.47 20.12 19.5 


从 表 3 可 以 看 出 ， 修 正 前 仿真 值 与 实验 值 的 最 
大 误差 为 20.12 "C. 修正 目标 函数 如 式 (3) 所 示 。 修 
正 过 程 如 下 , 首先 ,以 43、 互 、Bs 为 修正 参数 , 采用 
多 岛 群 遗传 算法 进行 修正 , 修正 前 误差 为 20.12 so 


修正 后 误差 为 9.13 °C, 2 ep 得 到 的 参 关 2 
作为 初始 值 ， 同 时 以 BB 人 


法 对 前 两 步 所 得 到 的 7 个 参数 进行 修正 ,修正 后 误 
差 由 6.67 °C 降低 到 了 1.14 "C。 各 参数 修正 前 与 修 
正 后 的 结果 如 表 4 所 示 。 

表 4 是 对 工作 在 3.6 A 电流 时 的 热 模 型 进行 的 
修正 ， 采 用 同样 的 修正 方法 ， 分 别 修正 工作 电流 为 
0.6 A、1.2 A、1.8 A、2.4 A、3.0 A 时 的 热 模型 .得 到 
修正 后 的 输入 参数 -电流 的 变化 曲线 , 如 图 5 所 示 。 

5(a) 中 ，N 性 材料 的 导热 系数 的 取 值 范 
围 由 原来 的 0.1~1.5 W/(m.*C) 变 为 0.95~1.45 
W/(m*C)。 赛 贝克 系数 取 值 范围 由 原来 的 -4 x 
10-4 ~ 一 1.2 x 10-4 V/°C 变 为 -2.07 x 10-4 ~ 
一 1.75 x 10-4 V/"C。P 型 材料 的 导热 系数 与 N 型 
相同 ， 赛 贝克 系数 随 电 流 增 大 而 减 小 ， 修 正 后 变化 
范围 为 10-4 ~ 1.25 x 10-4 VJ/*C。 导 热 硅胶 

rage 围 从 修正 之 前 前 的 1~3 W/m.*C 变 

怪 EE 的 1.65~2.85 W/(m*C)。 修正 后 的 热 端 对 

热 系数 取 值 随 着 电流 增 大 而 减 小 ， 取 值 范围 为 
25~54 W/(m2.°0), 冷 端 对 流 换 热 系 数 随 着 电流 增 大 
而 减 小 ， 由 修正 前 的 5~10 W/(m?.*C) 变 为 修正 后 
的 8.5~9.8 W/(m?.*C)。 参数 取 值 修正 前 后 发 生 了 明 
显 的 变化 。 


人 表 4 参数 修正 结果 
区 Table 4 Parameter correction result 


参数 修正 前 第 一 步 二 步 第 三 步 
A2/(W/(m-:°C)) 1.45 1.45 X_。 1.47 1.05 
A3/(V/°C) 一 1.38 x 10-4 -2.32x10-4 。 @~ _2.32 x 10-4 一 1.98 x 10-4 
B2/(W/(m:°C)) 0.6754 0.6754 1 1.22 
Bs/(V/°C) 0.72x10-4 0.79x10 A 0.79x10-4 0.79x10-4 
E/(W/(m-:°C)) 1.56 RC 1.38 7 
G1/(W/(m? .°0)) 53.2 Ce 38.4 25.2 
G2/(W/(m? .°C)) 5 SA。\5 5.36 8.75 
Object/°C 21 xw 9.13 6.67 上 及 
一 1.6 


N 导 热 系数 /W.m-LK-1 


05 10 15 20 25 30 35 40 
电流 /A 
(a) N 导 热 系数 随 电流 的 变化 曲线 


N 赛 贝克 系数 /10-4V.K- 
| | | | 
So Ne) co | 


| 
Dd 
~ 


| 
M2 
b 


05 10 15 20 25 30 35 4.0 
电流 /A 
(b) N 赛 贝克 系数 随 电流 的 变化 曲线 
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P 导 热 系数 /W.m-LK-! 


0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 
电流 /A 
(oO)P 了 导热 系数 随 电流 的 变化 曲线 


导热 硅胶 片 导热 系数 /Wim-K- 


05 10 1.5 2.0 N30 3.5 40 


由 
@ ee 


10.0 


a 
2o 


系数 /W.m2.K-1 
ed 
小 


9.2 


冷 端 对 流 换 热 
% 
Oo 


03 ”10 .18 


P 赛 贝克 系数 /10-4V.K-1 


05 10 15 20 25 30 35 4.0 
电流 /A 
(d) P 赛 贝克 系数 随 电流 的 变化 曲线 


05 10 15 20 25 30 35 4.0 
电流 /A 
( 热 端 对 流 换 热 系数 随 电流 变化 曲线 


2.0 2.5 3.0 3.5 4.0 


电流 /A 


(g) 冷 端 对 流 换 热 系数 随 电 流 变 化 曲线 
图 5 0.6~3.6 A 电流 下 各 参数 随 电 流 变 化 曲线 


Fig. 5 Temperature contrast of each point from 0.6~3.6 A 


由 于 将 所 有 影响 模型 仿真 精度 的 因素 简化 在 以 
上 7 个 参数 中 ， 且 这 7 个 参数 并 不 独立 ， 相 互 影 
响 。 使 得 修正 后 得 到 的 参数 曲线 不 再 具有 单一 的 热 
物性 参数 ， 所 以 得 到 的 参数 修正 曲线 变化 具有 明显 
的 非 线性 特征 .修正 得 到 的 参数 曲线 仅 作为 当前 仿 
真 模型 的 准确 输入 ， 据 此 建立 热电 制冷 器 的 仿真 模 


型 并 验证 其 模型 准确 性 . 取 电 流 值 分 别 为 0.9 A、1.5 
A、2.1 A、2.7 A、3.3 A 时 ,各 输入 参数 在 修正 后 对 
应 的 数值 , 结果 如 表 5 所 示 ， 

将 参数 带 入 有 限 元 模型 中 , 计算 出 a 点 .5 点 、c 
点 、d 点 的 温度 , 并 进行 实验 验证 。 实验 与 仿真 结果 
对 比如 图 6 所 示 ， 
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表 5 修正 后 参数 的 取 值 


Table 5 The revised parameter values 


电流 /A 0.9 1.5 
A2/(W/ (mm.°©0)) 1.3 1.0 
A3/(V/°C) 一 1.75BE-4 一 1.7E-4 
B2/(W/(m-:°C)) 1.32 1.4 
Bs/(V/°C) 1.2E-4 1.2E-4 
E/(W/(m:°C)) 2.48 2.25 
G1/(W/(m?:°C)) 49 41.5 
G2/(W/(m?:°C)) 9.4 9.05 


从 图 6(a) 可 以 看 出 , 测 点 a 的 实验 值 随 着 工作 
电流 的 增 大 而 增 大 , 仿真 值 随 着 电流 的 增 大 而 增加 ， 
且 趋 势 基本 一 致 . 修正 前 仿真 值 整体 小 于 实验 值 , 电 
流 为 0.5 A 时 , 误差 最 小 , 随 着 电流 增 大 而 增 大 , 电 


45 
一 口 一 修正 前 仿真 值 
40『| 一 Xx 一 修正 后 仿真 值 
一 A 一 实验 值 
35 
起 30 Oo 
us| 
25 :OD 
E 4 
20 ml 
15 
05 10 15 20 2.5 


3.0 
电流 /A 。 和 XT 
(a) a 点 修正 前 、 修 正 后 仿真 计算 温度 实测 温度 值 对 比 


100 


| 一 o 一 修正 前 仿 
80| 一 x 一 修正 后 仿真 


一 人 一 实验 值 
60 上 
荐 
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(b) 5b 点 修正 前 、 修 正 后 仿真 计算 温度 值 与 实测 温度 值 对 比 


40 


| 一 一 修正 前 仿真 什 
30| 一 x 一 修正 后 仿真 值 
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I 20r 
ee 
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EE a 
由 10F 
i 
| 人 本 
ot XXX 一 
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(b) c 点 修正 前 、 修 正 后 仿真 计算 温度 值 与 实测 温度 值 对 比 


2.1 2.7 3.3 
1.15 1.45 1.25 
一 1.85E-4 一 1.8E 一 4 —1.8E—4 
1.29 0.9 1.0 
1.0E-4 9E-5 9.1E-5 
2.55 2.2 1.8 
33 27 26.5 
9.4 9.0 8.5 

40 


一 5 一 修正 前 仿真 


30 上 | 一 x 一 修正 后 仿真 
| | 一 A 一 实验 值 严 
20 上 


Pi 


温度 /CC 
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、\ 


2 A 
稚 p=” 
CN x 


SS 05 10 15 20 25 3.0 35 40 
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(d) q 点 修正 前 、 修 正 后 仿真 计算 温度 值 与 实测 温度 值 对 比 
图 6 仿真 值 与 实验 值 的 对 比 


Fig. 6 Comparison of measuring value and 


experimental value 


流 为 3.9 A 时 , 误差 最 大 值 为 9.9 °C. 修正 后 仿真 什 
与 实验 值 的 误差 整体 减 小 , 在 3.9 A 处 最 大 值 为 4.9 
°C, 从 图 6(b) 可 以 看 出 , 测 点 5 温度 随 着 工作 电流 
的 增 大 而 增 友 , 念 货 值 随 着 电流 的 增 大 而 增 大 , 趋势 
一 致 . 电流 外 于 2.5 A 时 , 修正 前 的 仿真 值 与 实验 值 
误 益 镁 和, 大 于 2.5 A 时 ,误差 较 大 , 且 随 着 电流 增 
大 , 电流 为 3.9 A 时 误差 达到 13 *C, 修正 之 
,电流 小 于 3.5 A 时 , 仿真 值 与 实验 值 误差 小 于 1 
°C,，3.9 A 时 , 误差 为 2.5"C.。 从 图 6(c) 发 现 , 电流 
从 0 A 到 4A 变化 时 , 测 点 < 的 实验 值 与 仿真 值 出 
现 先 减 小 后 增 大 的 情况 .修正 前 仿真 值 与 实验 值 的 
误差 整体 较 大 , 电流 从 0.6 A 变化 到 3.9 A 时 , 误差 
从 4.6 °C 增 大 至 20.9 °C. 修正 后 仿真 值 与 实验 值 的 
误差 虽 现 出 中 间 大 两 边 小 的 趋势 ， 但 最 大 误差 仅 为 
2.2 °C. 从 图 6(d) 可 知 , 测 点 4 的 实验 值 和 仿真 值 变 
化 趋势 与 测 点 C 一 致 . 修正 之 前 误差 随 着 电流 增 大 
而 增 大 , 在 3.9 A 时 最 大 误差 为 20.6 °C, 修正 后 变 
为 3220。 
通过 对 比 4 个 测 点 的 仿真 值 与 实验 值 误差, 发 
现 最 大 误差 不 超过 4.9。C. 在 工作 电流 为 0~4 A 时 ， 
有 限 元 模型 仿真 值 与 实验 值 的 误差 波动 小 ， 稳 定性 
好 。 说明 修 正 后 的 输入 参数 准确 度 较 高 ， 进 而 有 效 
的 证 明了 热电 制冷 器 建 模 精 度 较 高 ， 
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3 结 论 

针对 热电 制冷 器 热 仿真 模型 精度 低 的 问题 ， 本 
文 开展 了 模型 建立 与 修正 技术 的 研究 工作 ,结果 显 
示 : 

1) 在 运用 Spearman 相关 级 数 对 参数 进行 灵敏 
度 辨 识 的 基础 上 , 分 析 得 出 N 型 材料 和 了 型 材料 赛 
贝克 系数 、 导 热 硅胶 片 导热 系数 为 模型 的 强 相 关 参 
数 ， 而 对 流 换 热 系数 、N 型 材料 和 P 型 材料 导热 系 
数 为 模型 的 弱 相 关 参 数 。 其 他 参数 相关 度 较 低 ， 可 
不 作为 修正 参数 

2) 采用 全 局 优化 算法 和 直接 搜索 优化 算法 ,并 
结合 试验 数据 修正 热 模型 ， 能 够 有 效 提 高 模型 的 仿 
真 精度 ， 

3) 基于 修正 后 的 热 模型 ， 对 比分 析 了 热电 制冷 
器 在 工作 电流 为 0.3~3.9 A 范围 内 ,模型 的 仿真 值 
与 实验 值 之 间 的 误差 。 发 现 所 有 测 点 的 最 大 误差 不 
超过 4.9 °C, 较 修正 前 的 20.9。C 有 明显 降低 , 表明 
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